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QCD two-loop amplitudes for e^e 3 jets: the fermionic contribution 
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We review a new technique for the calculation of two-loop amplitudes and discuss as an example the fermionic contributions 
to —> qqg. 


1. Introduction 

The construction of fully differential next-to-next- 
to-leading order (NNLO) programs is a challeng¬ 
ing problem and urgently needed to match the accu¬ 
racy reached in today’s collider experiments. While 
for inclusive observables like for example the total 
hadronic cross section in annihilation the step 
from NLO to NNLO is far less demanding - and has 
been taken a long time ago - the contrary is true for 
less inclusive quantities as for example jet rates. Es¬ 
sentially one has to adress two problems in jet physics 
at NNLO accuracy. The hrst one is related to the fact 
that one needs to treat two-loop integrales with more 
than one scale. These integrals are much harder to 
solve than one-scale integrals. The second problem 
which needs to be solved is the combination of vir¬ 
tual and real corrections - which are separately in¬ 
frarot (IR) divergent - into a IR hnite cross section. 
Given these difficulties it is clear that only for specihc 
reactions NNLO calculations can be envisaged. Im¬ 
portant reactions for which one should go beyond the 
NLO approximation are for example; Bhabha scatter¬ 
ing, pp,pp 2 jets, and ^ 3 jets. Among these 
examples the annihilation into 3 jets is of partic¬ 
ular interest. Historically this reaction was one of the 
first clean tests of Quantum Chromodynamics (QCD) 
as the theory of strong interaction. Today the 3-jet 
production in annihilation is still an excellent 

laboratory for precise tests of QCD and jet physics. In 
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particular —s- 3 jets is also a perfect reaction for 

measureing the QCD coupling with high accuracy. 
It is exactly here where the need for NNLO predic¬ 
tions for 3-jet production becomes most prominent: 
at present the measurements of are plagued by 
scale uncertainties of the theoretical predictions [jl|]. 
These uncertainties arise from uncalculated higher or¬ 
der contributions. Although formally of higher order 
these uncertainties can still be large. They can be 
reduced only by a NNLO calculation. The residual 
scale dependence of the 3-jet rate f$ can be studied 
using the following formula 
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Here pOjPi denote the coefficients of the QCD p- 
function, y is the resolution of the jet algorithm and 
A, B, and C are the LO, NLO, NNLO coefficients at 
= s. The residual scale dependence is illustrated 
in hg. 0 for the Durham algorithm at y = 0.1. For 
the NNLO curves three different scenarios have been 
assumed; one scenario in which the NNLO correc¬ 
tions vanish for = s, and two scenarios where a 
± 5 % correction at p^ = s have been assumed. It 
is clearly visible that the NNLO corrections improves 
the residual scale dependence. For a 3-jet prediction 








2 





Figure 1. Residual scale dependence of fi{y, y)- 


at NNLO accuracy several ingredients are needed; 
First the matrix elements for ^ QQggg, qqq'q'g 
have to be evaluated. They are obtained from a 
leading-order calculation and are known for a long 
time Second the one-loop matrix elements for 

^ qqgg,qqq'^ are needed. They have been 
calculated in refs. Finally one needs the 

NNLO amplitudes for —*■ qqg [pl p[ [TTI, [l^ [l^] . 

In the following the calculation of fermionic contri¬ 
bution of the NNLO matrix elements is reported. One 
should keep in mind that the combination of these 3 
ingredients into an infrared finite cross section is still 
a highly non-trivial - and so far unsolved - problem. 

2. A new approach to calculate NNLO amplitudes 

Due to the tremendous activities in the field enor¬ 
mous progress has been made in the past. All the 
relevant double box scalar integrals have been calcu¬ 
lated [|^[^|^|^. For a few reactions these re¬ 
sults have been used to obtain NNLO scattering am¬ 
plitudes To¬ 

day one can say that a standard approach to perform 
these calculations exist. Essentially it consists of three 
steps: 

1. Generation of relevant Feynman diagrams. 

2. Reduction of tensor integrals to scalar integrals 


via Schwinger parametrization 

3. Reduction of scalar integrals with raised pow¬ 
ers and in higher dimensions to master integrals 
via partial integration and Lorentz invariance 
equations 

The final result is than obtained in terms of master in¬ 
tegrals. They have to be calculated analytically. The 
bottleneck of this approach is that one has to create 
and to solve a huge system of equations to obtain the 
desired reduction to master integrals. In contrast to 
the naive expectation, topologies which at first glance 
look not so difficult may involve much more work 
than the more complicated once. For example it is 
almost trivial to find the reduction scheme for the pla¬ 
nar double box although the master integral itself is 
one of the most complicated once. In ref. [ ^ ] we 
have proposed a different method. The basic idea is 
that one applies only obvious reductions like for ex¬ 
ample the triangle rule. In particular one does not 
perform a complete reduction to master integrals. In¬ 
stead one calculates the scalar integrals in higher di¬ 
mensions and with raised powers of the propagators 
directly. To illustrate the method let us consider the 



Figure 2. Tensor reduction of the general pentabox 
topology a) to scalar integrals in higher dimensions 
and raised powers of the propagators b). Application 
of the triangle rule yielding the C-topology C{m = 

|-fe,Vi,V2,V3,V4,V5) C). 














3 


so-called penta-box. A corresponding Feynman di¬ 
agram is shown in fig. Applying Schwinger 

parameterization one arrives at the scalar topology 
shown schematically in fig. An obvious appli¬ 
cation of the triangle rule allows the elimination of 
two propagators [ |^ ] (cf. fig. ^). The resulting 
topology is often called C-topology. As mentioned 
earlier the C-topology appears in higher dimensions 
{d = 2m — 2e, m S N) and with raised powers of 
the propagators (vi,...,V 5 ). In the approach advo¬ 
cated here an analytic expression for the C-topology is 
needed. Such an expression can be obtained in terms 
of infinite sums The generic structure is shown 
in eq. (|). 
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It is obvious that such a representation is completely 
useless unless one is able to evaluate the e-expansion 
in terms of known functions. Fortunately this is possi¬ 
ble using a few basic algorithms for nested sums [|^. 
The basic idea here is that one rewrites the sums using 
the relation 


harmonic sums by 


n ii ‘k-1 ^'i ^‘k 

S{n-,mi,...,mk-,xu--.,Xk) = ^ 

il = ll2 = l 
n 

= L -:^S[i\m2,...,mk\X2,...,Xk) (5) 

/=l ^ 


the following four basic operations are sufficient to 
evaluate the e-expansion of the C-topology in terms 
of multiple poly logarithms [ |T^ ] : 

1. Multiplication 


S{n\mi,...\xi,...) X S{n\m\,...\x\,...) (6) 


2. Sums involving i and n — i: 


n—l 

E ■^^S{i-m2...-,X2,...) 
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3. Conjugations: 

E (^ ” ^{-'^Y^S(i\mi,...,mk\x\,...,Xk) (8) 

4. Sums involving binomials and n — i: 
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Furthermore one can expand the F-functions by the 
means of 

r(-n-be) = ■tlEr(H-e)(l-b5i(n)e 

enl 

+5ri(n)e^ -l-5i,i,i(n)e^ -l-5i,i,i,i(n)e"* + • • • (4) 

and related identities [|^. Here ...(«) denote the 
harmonic sums, see e.g. ref. [|^]. Inspecting the sums 
which appear using eq. and eq. (Q) one observes 
that the algorithms needed are just a generalization 
of the algorithms studied in ref. [Q]. In particular 
defining the nested sums as a generalization of the 


/ n—i 

^ r W S{n-i-m2,...-,X2,-) (9) 

In ref. [ |^ ] we have given explicit algorithms for 
the four basic operations. For their implementa¬ 
tion we used Form3 [|^ and Ginac Us¬ 

ing these programs we calculated more than 300 
C-topologies analytically. In particular we have 
checked that our results for the two master in¬ 
tegrals C(2,1,1,1,1,1),C(2,1,1,1,1,2) agree with 
those given in ref. All the simpler topologies 

can be calculated in the same way. 
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3. Results 


The amplitude for the process ^ y* ^ 

Q{pi)s{P 2 )q{p^) can be written as a leptonic part 
multiplied by a hadronic part 





with s = {p\+ P 2 + Pi)^ and Further¬ 

more the hadronic part can be decomposed according 
to the colour structure. In particular we have 
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with the generator of the SU (N) gauge group and 
tif the number of massless quarks. In this talk we dis¬ 
cuss only the contributions proportional to tifN and 
The hadronic contribution can also be decom¬ 
posed according to the spinor structure: 

Hf, = cP{pi\^2\P3)£gf, 

+ ■ ■ ■ + C4^{Pl\h\P3){^g ■ Pl)Pl^, + ■ ■ ■ 

+Cl3^{pi\h\P3){^g-P3)P2^, ( 12 ) 

where the functions ci — C 13 depend only on the ra¬ 
tios xi = 2{p\ ■ P 2 )/s, X 2 = 2{p2 ■ P 3 )/s. Due to vari¬ 
ous constraints arising from gauge invariance or quark 
anti-quark antisymmetry only 4 of them are needed. 
All the remaining ones can be expressed in terms of 
these 4 functions which we chose to be C 2 , C 4 , eg and 
C 12 . The perturbative expansion in a., of the functions 
a is defined through 


Ci = g!i 1 ef'-f 


(|)/.4(t|)/f4<„„3)-<.3. 


Using the methods described in the previous section 
we calculated all the 13 functions and used the con¬ 
straints as a check of our calculation. Furthermore 
we checked the structure of the ultraviolet as well as 


the infrared divergencies Note that the gen¬ 

eral structure of the IR divergencies as conjectured 
in ref. has been proven recently by Sterman and 
Tejeda-Yeomans [pT|]. We have also compared our 
results with results published recently [|^||]. In refs. 

^ the result is written in terms of 2 -dimensional 
harmonic polylogarithms instead of the multiple poly¬ 
logarithms used here. Expressing the 2-dimensional 
harmonic polylogarithms in terms of multiple poly¬ 
logarithms we found complete agreement for the nfN 
and ^ contribution calculated here. 

Following ref. [^] we used the general structure of 
the infrared divergencies as described by Catani [ p^ 
to define the finite parts of the coefficients cf. 

^(2),fin ^ cf>’'“-l(')(e)4‘^"“-l(2)(e)4'’\(14) 

with the one- and two-loop insertion operators l(^) (e) 
and l( 2 )(e) given in ref. [p^. As an example, we 
present the result for the n/A-contribution to the fi¬ 
nite part at two loops, 
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We have introduced the function R(jii:i ,;iC 2 ) defined in 
[p^. In addition, it is convenient, to define the sym¬ 
metric function Ri (xi ,^ 2 ), which contains a particular 
combination of multiple polylogarithms [p5|]. 


Ri(xi,X 2 ) = ( ln(xi)LiLi 


XI 


,Xl-|-X2 


X1-I-X2 

-i^2ln(l-xi-X2) - iln(xi)ln(x2)ln(l-xi-X2) 


-ln(xi)Li 2 (xi-|-X 2 ) -Lip 2 

Xl 


XI 


-,X1-|-X2 


—Lh.i 


X1-I-X2 


X1-I-X2 

,Xi-|-X2 ) -|-Li3(xi-|-X2) ) + (xi ^X2). 


All the multiple polylogarithms have simple argu¬ 
ments. As a consequence it is straightforward to ob¬ 
tain the analytic continuation, which can be used for 
the crossing of the amplitude. 


4. Conclusion 

In this talk we have presented the calculation of a 
specific colour structure of the QCD two-loop ampli¬ 
tude ^ qqg. The computation has been done 
by a completely new method . The new approach 
provides a very efficient method for two-loop calcu¬ 
lations. The tools developed in this approach can also 
be used for a systematic expansion of generalized 
hypergeometric functions. In addition the presented 
results provide an important cross check on recently 
obtained results [^,^. 
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